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Abstract. We present a model in which AGN jets are described by a dense series of plasma blobs, distributed in 
masses and velocities. The blobs expand and collide, and ultimately form a continuous jet flow as observed on larger 
scales. When blobs collide, energetic electrons are produced by first order Fermi acceleration, and emit synchrotron 
radiation at the radio-optical frequencies with a double-broken power law spectrum. The difference of this model 
to the well-known Marscher&Gear-model is, that in our model shocks are produced on a range of scales, and exist 
only temporarily, which causes differences in the prediction for the spectral evolution of flares. We apply our model 
to radio-submm data obtained for 3C454.3, and briefly discuss implications for gamma-ray production. 



1. Introduction 



Based on the work of Marscher fc Gear (1985)| flares in 
AGN jets usually modelled by shocks propagating in a 
continuous plasma jet from its origin to the outside. While 
the shock continuously energizes electrons, synchrotron 
radiation emitted by it goes through a flare, owing to the 
competing cooling processes of inverse Compton radiation, 
dominating the in the rise phase, and adiabatic losses in 
the expanding jet, dominating in the decay phase, with 
the possibility of a plateau phase in-between, in which 
synchrotron radiation itself is the dominating cooling pro- 
cess. As the relative strength of cooling processes depends 
on global jet parameters, which are unlikely to change 
on a short timescale, all flares should evolve in a similar 
way, with their rise, plateau and decay timescales being of 
the same order of magnitude, typically weeks to months. 
While this fits observations of radio flares quite well, the 
situation is less clear in the optical and gamma-ray regime. 
Here, flares appear on a large range of time scales, from 
less than an hour in some TeV blazars, to weeks for strong 
blazars in the Fermi regime. 

To accommodate the situation in the gamma-ray 



regime, Spada et al. (2001) suggested a so-called "inter- 
nal shock model" for blazars, motivated by the model 



of Rees & Meszaros (1994) for the explanation of rapidly 
variable gamma-ray burst lightcurves. In this model, 
shocks are formed in random collisions of shells, and 
subsequently accelerate electrons to emit synchrotron 
and inverse Compton radiation. Within the frame- 
work of this synchrotron-self Compton (SSC) scenario, 



Spada et al. (2001) used their model to explain qualita- 



tively the global characteristics of blazar variability at op- 
tical to gamma-ray frequencies. 

In the current work we revoke the internal shock model 
for blazars, focussing on the explanation of flare evolution 
in the radio-submm regime, which we claim to allow us to 



pinpoint the physical parameters of the jet independent 
on the assumptions for gamma-ray production. A detailed 



discussion will be given by |Haberlein fc Rachen (2010) 



henceforth referred to as "Paper I" . Here we discuss briefly 
the basic idea and properties of this model, and apply it 
to a flare observed in the bright blazar 3C454.3 in 2005. 

2. Jets, shocks and plasmons : not a contradiction 

Essentially two fundamental mechanisms have been pro- 
posed to explain why plasma jets may be launched from 
accreting black holes: the [Blandford fc Znajek (1977) pro- 



cess, which is based on electromagnetic extraction of the 
angular momentum of a rotating black hole in form of a 
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strong Poynting flux, and the Blandford & Payne (1982) 
process of hydromagnetic extraction of rotating plasma 
from the corona of the accretion disk. In the recent 
discussion of AGN and other accreting sources, it is 
generally assumed that both mechanisms contribute to 
the jet formation, possibly with different weights de- 
pending on the physical properties of the accreting sys- 
tem [e.g., Meier et al. (2001)] )]. Although plasma physi- 
cists have shown that solutions for the continuous emission 
of jets exist, all involved physical processes are prone to in- 
stabilities, in particular in the highly relativistic regime of 
accreting black holes. This justifies the assumption that, 
although jets are observed as continuous plasma flows on 
larger scales, they are not born like this, but as a dense, 
random flow of plasma blobs, emitted a central source 
similarly to a discharge process. 

Let us assume that the blobs are emitted on an aver- 
age rate /* < c/Rs, if Rs is the Schwarzschild radius of 
the central black hole, at highly relativistic speeds with a 
bulk Lorentz factor T* ~ 10. Being generally unconfined, 
in their comoving frame the blobs would expand in all 
directions approximately with their speed of sound, /3 s c, 
leading to an opening angle of the yet discontinuous jet 
of /3 S T~ . In the initial phase, where relativistic gas domi- 
nates the pressure, the expansion decreases with distance 
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Fig. 1. Illustration of the discontinuous jet model. Plasmons created some 10—100 Schwarzschild-radii from the black hole 
collide mostly at around 10 4 -Rs, converting a significant fraction of the bulk Lorentz factor variance into internal energy, and 
ultimately into radiation. As the "collision region" is quasi-stationary in position and emits most of the radiation in a spectrum 
averaging out the individual collisions, it might be identified with the bright "millimetre-core" seen in many blazars. Particularly 
large differences in bulk Lorentz factor may lead to early, strong collisions which appear as flares in the spectrum. As the collision 
region is leveling out the differences in bulk speed, the jet merges afterwards quite smoothly. A series of fast blobs, e.g. after a 
surge in the accretion rate, may propagate through the whole system, leading to a series of flares, and ultimately form a moving 



shock in the continuous jet, explaining VLBI observations as reported by Marscher et al. (2008) 



r from the black hole approximately as A oc r -1 / 2 , leading 
to an apparent collimation of the jet to angles of less than 
1°. After this, we assume conical expansion with A ~ 0.1. 
Neglecting velocity differences of the blobs, they would 
merge rather smoothly after a distance r mrg ~ r 2 c//*A, 
as their expansion has filled up the gaps between them. 
Of course, a realistic "discharge" process near the black 
hole will emit the blobs with randomly distributed Lorentz 
factors + 5Ti, with (STi) = and a normalized stan- 
dard deviation or = \J /T*. Inelastic collisions be- 
tween plasmons will then become frequent at a distance 
^coi ~ r*c//*crr, converting on average a fraction Up/4 of 
the bulk flow into internal energy. Subsequent collisions 
will lead to more and more elongated "jet-bits" , with the 
spread in Lorentz factor diminishing as ar' fc ' ~ ar^/2 k , 
if k denotes the number of collisions an individual blob 
has suffered. Eventually, <C As, and blob-expansion 

is again the dominating process to form the continuous jet 
at > 10 5 i?s- Figure [l] illustrates this process. 

We want to emphasize that our plasmons cannot be 
identified with the observable VLBI components. In our 
model, there are ~ r*/A ~ 10 3 plasmons in the discon- 
tinuous phase of the jet (r < 10 5 i?s) at any time. They 
are usually dark, as they are born with extremely high 
magnetic field cooling down any synchrotron emission on 
a short timescale. The blobs become bright only after they 
are re-energized at ~ 10 4 — 10 5 i?s m collisions. For most 
blazars, even the best VLBI pictures will not be able to 
resolve single collisions; rather we expect that the known 
resolved VLBI features in AGN jets correspond to the dis- 
tinct zones in our model, as shown and explained in Fig.[l] 
We also emphasize that, in general, collisions in our model 
do not correspond to a "flares" in the observed spectrum. 
As the rate of collisions must be equal to /* to create 
a stationary system, that means for /* ~ c/10i?s and a 
10 8 M Q black hole several collisions per day. As the syn- 



chrotron cooling time scale is much larger, i.e. f*t syn » 1, 
we see at any instant the superposed emission of many 
plasmon collisions, averaging out their spectral and tem- 
poral properties. The average collisions with ST /T < or 
therefore represent what observes call the quiescent emis- 
sion of a blazar; flares, in the contrary, are collisions with 
8T/T > 3fJr, which can have luminosities up to three or- 
ders of magnitude higher than the average, as we will show 
below. In Paper I we present a full scale simulation of the 
jet, starting from the statistical properties of the central 
engine and calculating of all n-th order collisions on the 
way to the continuous jet. There we also we demonstrate 
that the overlay of average "mini-flares" indeed produces 
a broken power law spectrum with a sa below, and 
a m 0.7—1.2 above the break frequency typically in the 
range 10— 300 GHz, as usually observed in blazars. 

3. Flare characteristics in a discontinuous jet 

When two plasma blobs, both carrying an energy ss Mc 2 F, 
collide with a Lorentz factor difference 5T, two shocks run- 
ning in opposite directions with speeds Ah ~ \/2STJT 
(for 5T <C r) are formed, and internal energy E- lnt ~ 
jMc 2 (ST /r) 2 is created. We assume that a significant 
part of this internal energy is going into turbulence, specif- 
ically magnetic field creation in turbulent dynamos, and 
another significant part into nonthermal particles through 
the mechanism of first order Fermi acceleration at the 
shocks. First order Fermi acceleration has the property 
to produce power law particle distributions, and there- 
fore to explain the ubiquitous presence of power law spec- 
tra in sources known to contain shocks (see Drury (1983) 
and references therein). We assume in the following the 
canonical power law index for Fermi-accelerated electrons, 
i.e. iV(7) = N e j~ s with s = 2, for the electron Lorentz 
factor in the range 7 < 7 < 7. We introduce the pa- 
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rameters £g and £ e by ug = B 2 /8n = ^gE- mt /2V and 



N e m e c 2 ln( 7 / 7 ) = £ e E int /2V, where V 



the volume of a "virgin" plasmon. In a conical jet, R re- 
lated to the distance r from the black hole as R « rf3 s /T. 
As the synchrotron luminosity L syn oc R 2 N e B 2 , we get 
L syn oc E 2 nt /R 4 oc (<5r/r) 8 for constant 7/7 [] Therefore, 
"3cr-events" in the distribution of collisions can outshine 
the sum of all average collisions, revealing temporarily the 
spectral properties of a single emission region with a spe- 
cific size i?, and specific values for plasma properties like 
N e or B. This is, what we usually call a flare. 

The spectrum of a flare is represented as a broken 
power law with three characteristic frequencies. The first 
is the synchrotron-self absorption frequency v a , defined by 
requiring that the opacity r ssa (t' a ) = 1, yielding 



<7 T RugU e 



Am e c 2 01(7/7) 



(1) 



where er T = 8irr 2 /3 is the Thomson cross section. For 
v < v a , the spectrum is S v oc j/ 5 / 2 , for v > v & we will have 
S v oc v~ a with a = (s — l)/2 = 0.5. This spectral index, 
however, is only correct for t syn (v) > t, i.e., synchrotron 
cooling is not yet efficient for the electrons emitting at a 
characteristic frequency 1/(7) = 3ei?7 2 /16m e c at a time 
t after the collision. Above a &h(t) for which t syn = t, 
the electron spectrum steepens as s — > s + I, yielding 
S v oc v^ 1 . Evaluating v\,{f) for the fully developed flare, 
i.e., t ~ R/cfish, we obtain 



Rc? 2 h 



eBR 



(2) 



Finally, the spectrum reaches a maximum frequency at 
v c = ^(7)1 above which it turns into an exponential 
cutoff. The value of 7 is hereby determined by the re- 
lation i acc (7) = t syn (j), where i acc (7) « 3/t/c 2 ^ is 
the acceleration time scale for first order Fermi acceler- 
ation. |Biermann fc Strittmatter (1987) connected the dif- 
fusion coefficient k to the plasma turbulence spectrum, 
and showed that for a Kolmogoroff spectrum known from 
fully developed hydrodynamical turbulence, v c is expected 
in the infrared to optical regime, as observed for many 
AGN. In Paper I, we will discuss a slightly modified ver- 
sion of the BS87 approach in more detail; here, we focus on 
the break frequencies v & and v^,, which are more relevant 
in the radio-submm regime. 

We evaluate v a and numerically for the "average" 
collision, i.e., with R = TwA^r, = T*c/3 S / f*a r , /3 sh = 
y/2a^, for values /3 S = 0.1, /* = c/I00i? s , and 01(7/7) = 
7.5. If we further assume that the jet power is equal to 
the Eddington luminosity of the black hole, with mass 
10 9 M 9 M„, we find 



"0> 

8I0 13 Hz 4tg /3 e e /3 Ms 1/3 --, 
1.5 I0 6 Hz a- 9/ % 3/2 M~ 1/2 



200 GHz/200 GHz (3) 
-> I50GHz/lGHz (4) 



1 |Spada et al. (2001)| vary 7 rather than N e to account for 
constant fraction of nonthermal energy in electrons to total 
internal energy, which would lead to L syn oc (<5IyT) . 



We defined 1/ = Dv to consider that blazar emission is 
usually Doppler boosted, setting D = r* = 10. The values 
after the arrow correspond to ap — 0.3 and M9 = 5 (as 
expected for 3C454.3), the first value for £g — 0.3 and 
£ e = 0.01, and the second for £g = 0.01 and £ e = 0.3 
(see Section [5] for a discussion). We see that for all cases 
Vh < v ai that means we expect in flares a steep spectrum 
oc v^ 1 above v a . Note, however, that during the rise phase, 
v(t) > z^b, so that the flare may temporarily exhibit a 



; -0.5 



spectrum directly above v & 



Of particular interest is the behavior of the flux at 
v a during the flare evolution, as this corresponds to the 
(t'm, S m ) plots often used in radio astronomy to describe 
flare evolution. In our model, B is approximately constant 
while the shock is active, and R hardly changes during the 
shock crossing time. Therefore, we expect that S m rises for 
constant v m . When the flare decays adiabatically after the 
shock ceases, the flux in the optically thick region evolves 
as S v (x R 2 B~ X / 2 , as the source function does not de- 
pend on the electron spectrum. With S u oc v°l 2 and for 



B oc R , one obtains S„ 



(one can easily check 



that S m oc is obtained for B oc R^ 1 ). If synchrotron 
cooling dominates over adiabatic cooling at v & , which is 
normally the case initially as <C i^ a , one can show that 

5 /2 

S m oc v>nl , i.e., the flare decays along the spectral edge 
of self-absorption. This may be compared with the pre- 



dictions of the Marscher & Gear (1985) model, which are 
quite different. A more detailed discussion of flare evolu- 
tion in our model will be given in Paper I. 

4. Fitting the 2005 flare in 3C454.3 

To illustrate the capabilities of our approach to reproduce 
the temporal-spectral evolution of blazar flares, we apply 
our model to data obtained on the bright blazar 3C454.3 in 
a campaign which has been started in 2005 after an initial 
strong mm-/optical flare. The main goal of the campaign 
has been to obtain quasi-simultaneous radio-spectra cov- 
ering the frequency range from I - 300 GHz, which can be 
used to study the spectral evolution over time. 

The Effelsberg 100 m radio-telescope of the Max- 
Planck-Institut fur Radioastronomie with its sequence of 
secondary focus heterordyne receivers is well suited to 
measure flux densities quasi-simultaneously and within 
30-40 mins at 2.7, 5, 8.4, 10.5, 15, 22, 32, and 43 GHz. 
These measurements were occasionally complemented 
by flux density measurements at 0.6-1.4 GHz, when 
the corresponding prime-focus receivers were available. 
The data were analyzed and calibrated in the stan- 
dard way following similar procedures as described in 
|Fuhrmann et al. (2008)] and |Angelakis et al. (2008)| The 
Effelsberg radio data were complemented by measure- 
ments within a few days of the millimeter flux densities 
obtained through the regular AGN monitoring program, 
which is performed at the IRAM 30 m telescope on Pico 
Veleta (Spain). This instrument provided measurements 



at 90 GHz and 230 GHz (see Ungerechts et al. (1998) and 



Agudo et al. (2006) for more details and a description 
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frequency (GHz) 

Fig. 2. Fit of a single flare in the discontinuous jet 
model to data of the 2005/06 flare from 3C454.3, see also 
Krichbaum et al. (2006)| and the discussion in the text. 



of the data reduction) . We further included 230 GHz 
and 345 GHz data from the AGN-monitoring of the sub- 
millimeter array (SMA) on Mauna Kea (Hawaii, see 
Gurwell et al. (2007) and refs. therein), which are pub- 
lished in part by Villata et al. (2009) and refs. therein. 

Figure [2] shows the data together with spectra simu- 
lated from our full model, as described in Paper I. We as- 
sumed that two plasmons with Mc 2 = 7 10 54 erg, r x = 10 
and T2 = 13 collide at a distance ~10 19 cm from the 
black hole, consistent with a jet power of 10 48 erg/sec 
and /* ~ 10 -6 Hz. We further assumed (3 S = 0.1, thus 
R r-j 10 17 cm, and a constant viewing angle of 5° to 
calculate the Doppler factor. All parameters have then 
been optimized in an adaptive method by minimisation 
of J2tJ2iXit f° r an data sets taken at the times t as 
given in the plot, where the simulated flare spectra have 
been added to a constant underground fitted to the data 
from September 2004. As the starting time of the flare is 
unknown, we have identified the data set from July 4th 
with the maximum of the flare. Frequencies and times 
have been properly corrected for the redshift of the source, 
z = 0.859. As the best fit parameters, we found £b ~ 0.3, 
and £ e « 10 -3 . The temporal evolution of the flare is rep- 
resented quite well, until presumably a new flare arises 
beginning of February 2006. The discrepancies at the low 
frequency turnover from the background to the flare spec- 
trum should not surprise, as our model expects the back- 
ground to change with time. Clearly, data at frequencies 
above 350 GHz would have been helpful to further con- 
strain our parameters. 

5. Implications for gamma-ray emission 



In their discontinuous jet model, Spada et al. (2001) in- 
cluded inverse Compton (both SSC and external IC) emis- 
sion, and determine the necessary jet parameters from the 
requirement to fit existing gamma-ray data. In contrast, it 
is our approach to determine the jet parameters from syn- 
chrotron emission alone, and keep minds open regarding 



the mechanism of gamma-ray production. We may remind 
the readership at this point that there are alternatives 
to inverse Compton models for gamma-ray production. 
In the so-called "hadronic" scenario, protons are acceler- 
ated to extremely high energies, enabling them to produce 
gamma-rays by pion induced cascades Mannhei m (1993)] , 
with a possible admixture of proton-synchrotron radia- 
tion Miicke et al. (2003)] . As first order Fermi accelera- 
tion, which is the base of our model, does not distinguish 
between protons and electrons in principle, we may state 
that hadronic models are at least as "natural" as the "lep- 
tonic" inverse- Compton scenario. However, hadronic mod- 
els are very difficult to handle technically, so we abstained 
from making detailed predictions for the gamma-ray spec- 
tra in our model at first; nevertheless it is possible to draw 
some preliminary conclusions on gamma-ray production 
from our results. 



As discussed by Rachen (2000) leptonic and hadronic 
models can be distinguished by the required relation of 
magnetic field to photon density, i.e., by the ratio Cs/^e: 
while leptonic models need £s < £ e (in case of TeV blazars 
even £3 <C £ e ) to explain gamma-ray data, hadronic mod- 
els require £g ~ £p 3> £e, where £ p denotes the fraction 
of internal energy in relativistic protons. For our fit of the 
flare of 3C454.3, £b 3> £ e was found. It would be prema- 
ture to consider this evidence for hadronic origin of gamma 
rays from 3C454.3, as our fit method does not really per- 
form a full parameter space search, identifying all allowed 
and excluded regions. However, it does give confidence in 
our method to distinguish between hadronic and leptonic 
models in future, improved applications. This certainly re- 
quires to include all gamma-ray production methods into 
our model, and then compare with both Fermi gamma-ray, 
and detailed radio-optical data. For the latter, as obvious 
from the discussion in Section [4] closing the infrared gap 
in regular blazar monitoring would be highly desirable. 
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